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Abstract  systems  for  rotorcraft  transmissions  often  fail  to  provide 

sufficient  time  between  warning  and  failure  so  that  safety 
A  comprehensive  procedure  in  predicting  faults  in  gear  procedures  can  be  implemented.  In  addition,  inaccurate 

transmission  systems  under  normal  operating  conditions  is  interpretation  of  operating  conditions  often  result  in  false 

presented.  Experimental  data  was  obtained  from  a  spiral  alarms  and  unnecessary  repairs  and  downtoe.  The  develop- 

bevel  gear  fatigue  test  rig  at  NASA  Lewis  Research  Center,  ment  of  an  accurate  machine  health  monitoring  system  for 

Time  synchronous  averaged  vibration  data  was  recorded  rotorcraft  transmission  systems  will  increase  safety  and 

throughout  the  test  as  the  fault  progressed  from  a  small  reduce  ftnancial  investments, 

single  pit  to  severe  pitting  over  several  teeth,  and  finally 

tooth  fracture.  A  numerical  procedure  based  on  the  Wigner-  Three  major  types  of  vibration  analysis  have  been  used 

Ville  distribution  was  used  to  examine  the  time  averaged  in  the  early  detection  of  machine  failure,  namely,  frequency 

vibration  data.  Results  from  the  Wigner-ViUe  procedure  are  analysis,  various  techniques  using  the  time  domain  analysis, 

compared  to  results  from  a  variety  of  signal  analysis  and  joint  time-fisquency  analysis, 

techniques  which  include  time  domain  analysis  methods  and 

frequency  analysis  methods.  Using  photographs  of  the  gear  For  frequency  analysis,  the  ampUtude  spectrum  of  the 

tooth  at  various  stages  of  damage,  the  limitations  and  measured  vibration  signal  is  calculated  and  displayed  m  a 

accuracy  of  the  various  techniques  are  compared  and  continuous  manner.^^’^  This  method  is  particularly  power- 

discussed.  Conclusions  are  drawn  from  the  comparison  of  fill  in  identifying  vibrations  due  to  various  mechanical 

the  different  approaches  as  well  as  the  applicability  of  the  elements  in  the  system  at  their  operating  frequencies.  For 

Wigner-Ville  method  in  predicting  gear  faults.  example,  a  gearbox  will  experience  vibration  at  the  funda¬ 

mental  and  multiple  harmonic  frequencies  of  the  gear  tooth 
I.  Introduction  meshing  frequency,  which  are  usually  different  for  each  pair 

of  gears  in  the  gearbox.  The  changes  in  the  frequency 
In  the  last  two  decades,  with  the  ever  increasing  spectram,  such  as  an  increase  in  the  vibration  amplitude  at 
for  higher  power  and  performance,  premature  a  particular  frequency,  or  the  appearance  of  modulation 
failures  in  high  performance  machinery  often  result  in  major  sidebands  about  that  frequency,  can  often  be  identified  and 
production  and  finanrial  losses,  and  sometimes  even  lead  to  related  to  a  single  element  in  the  complex  system.  Further- 
catastrophic  consequences.  In  the  aerospace  industry,  one  of  more,  abnormal  vibration  signals,  caused  by  we^,  deforma- 
the  major  concerns  is  the  rotorcraft  transmission  system,  tion  of  the  machine  element,  "ghost  component  etc.,  which 
Currently,  on-board  operational  condition  monitoring  are  sometimes  difficult  to  quantify  in  the  time  domain,  can 
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be  readily  identified  in  the  fi-equency  spectrum.  However, 
the  spectral  analysis  technique  is  difficult  to  apply  in  a 
highly  con^lex  system,  such  as  a  helicopter  gearbox,  which 
may  have  as  many  as  30  gears  and  50  bearings.  The  large 
number  of  spectral  lines  which  these  elements  produce 
during  normal  operation  can  make  it  difficult  to  even  detect 
changes  in  the  spectram,  let  alone  to  identify  the  sources  by 
which  these  changes  are  generated.  In  addition,  the  modula¬ 
tion  of  the  vibration  signal  of  early  fatigue  cracks  in  the 
mechanical  elements  often  yield  sidebands  which  extended 
over  a  broad  frequency  range  with  low  amplitude.  Under 
these  circumstances,  the  background  noise  level  in  the 
spectmm  can  easily  mask  the  presence  of  the  sidebands,  and 
mechanical  defects  in  the  system  may  go  undetected. 

Time  domain  analysis  is  also  used  for  the  early 
detection  of  machine  faults.^*  The  digitized  time 

domain  vibration  signal  is  averaged  over  a  large  number  of 
cycles,  synchronous  with  the  ranning  speed  of  a  particular 
machine  element.  This  time  synchronous  averaging  removes 
the  background  noise,  leaving  a  more  accurate  estimate  of 
the  vibration  signal  of  the  component  of  interest.  This 
process  can  be  repeated  for  each  machine  element  in  the 
system  to  provide  the  vibration  signal  relative  to  that 
particular  element.  A  variety  of  statistical  based  techniques 
are  then  applied  to  the  averaged  signal  to  extract  specific 
fault  information.  In  most  machinery,  vibration  signals  are 
often  characterized  by  a  variation  of  both  amplitude  and 
frequency  with  respect  to  time.  For  most  mechanical  faults, 
both  the  amplitude  and  frequency  contents  of  the  signal  will 
shift  during  the  period  in  which  the  damage  is  generated. 
The  nature  of  these  short-term  variations  in  the  time  signal 
could  prevent  an  accurate  indication  of  the  characteristics  of 
the  mechanical  fault  present. 

The  joint  time-frequency  analysis  approach  applies  the 
Wigner-Ville  Distribution  (WVD)  1“3,14.16 
vibration  signal  of  the  system.  Unlike  the  Fourier  Transform 
process,  the  WVD  provides  a  instantaneous  frequency 
spectmm  of  the  system.  The  spectral  density  of  the  funda¬ 
mental  exciting  frequency  and  its  sidebands  changes  as  the 
shaft  rotates  through  a  complete  revolution.  If  the  vibration 
repeats  itself  each  revolution,  a  Fourier  Transform  will  show 
a  constant  frequency  spectmm,  while  the  changing  spectral 
density  with  time  from  the  WVD  provides  information 
concerning  the  frequency  distribution  concentrated  at  that 
instant  around  the  instantaneous  frequency.  Some  success 
has  been  achieved  in  applying  the  AWD  in  gear  transmis- 
Sion  systems  ’’’’  to  recognize  faults  at  various  locations 
of  the  gear. 


The  time  domain  analysis  includes  the  time  synchronous 
averaged  signal,  and  the  statistical  based  techniques  FMO, 
FM4,  NA4*,  and  NB4*  applied  to  the  time  averaged  signal. 
The  joint  time-frequency  analysis  uses  the  WVD  on  the 
vibration  data  with  special  windowing  techniques.  All  of  the 
analysis  methods  are  rqiplied  to  the  vibration  data  of  a 
failure  from  the  spiral  bevel  gear  fatigue  test  rig  in  the 
NASA  Lewis  Research  Center.  Results  from  the  various 
methods  are  compared  and  general  conclusions  are  drawn 
from  the  results. 

n.  Technical  Approach 

As  discussed  in  the  previous  section,  three  major 
methodologies:  A)  the  frequency  domain  ^proach,  B)  the 
time  domain  approach,  and  C)  the  joint  time-frequency 
approach  are  used  in  this  study.  The  following  is  a  descrip¬ 
tion  of  the  three  methodologies: 

(A)  Frequency  Domain  Techniques 

The  frequency  spectmm  is  found  by  applying  a  discrete 
FFT  on  the  time  averaged  signal  x(t),  such  that  the  spectral 
components  are 


where 

X(k)  =  Tj;  x(t)  expj'j^^' 

x(t) 

=  time  domain  signal 

X(k) 

=  frequency  domain  signal 

T 

=  sampling  time  interval 

N 

=  number  of  data  points 

The  frequency  components  are  examined  in  the  frequency 
domain  and  compared  with  those  obtained  at  various  stages 
of  the  fault  development  in  the  spiral  bevel  pinion. 

(B)  Time  Domain  Techniques 

Four  different  time  domain  techniques  for  early 
detection  of  gear  tooth  damage  are  used  in  this  study  for 
evaluation  and  comparison.  All  of  the  time  domain  tech¬ 
niques  are  ^plied  to  the  vibration  signal  after  it  has  been 
time  synchronously  averaged.  These  techniques  are:  FMO, 
FM4,  NA4*,  and  NB4*.  These  parameters  are  defined  as 
follows: 

FMO*^  is  a  course  fault  detection  parameter  that 
conq}ares  the  maximum  peak  to  peak  amplitude  to  the  sum 
of  the  mesh  frequencies  and  its  harmonics 


The  objective  of  this  paper  is  to  examine  and  compare 
the  three  different  approaches  to  detect  gear  wear  and 
failure.  The  frequency  domain  analysis  is  based  on  the 
spectral  display  from  a  Fast  Fourier  Transform  algorithm. 


FMO  = 


PP 

E  A(fi) 

i=l 


(2) 
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where 

PP  =  maximum  peak  to  peak  amplitude  in  signal 
A(fj)  =  amplitude  at  mesh  frequency  and  harmonics 
n  =  totad  number  of  harmonics  in  frequency  range 


FM4‘'^  is  an  isolated  fault  detection  parameter,  and  is 
given  by  the  normalized  kurtosis,  of  the  resulting  difference 
signal  as 


FM4  = 


(di  -  d)'‘ 


(3) 


where 

d(t)  =  A(t)  -  R(t) 

d  =  mean  value  of  d(t) 

A(t)  =  original  time  synchronous  signal 
R(t)  =  regular  meshing  components  plus  their  first 
order  side  bands 

N  =  total  number  of  data  points  in  the  time  signal 


NA4*^’^°’^^  is  a  general  fault  detection  parameter  with 
trending  cjpabilities.  A  residual  signal  is  constructed  by 
removing  regular  meshing  components  from  the  tune 
averaged  signal.  For  NA4,  the  first  order  sidebands  stay  in 
the  residual  signal  and  the  fourth  statistical  moment  of  the 
residual  signal  is  then  divided  by  the  averaged  variance  of 
the  residual  signal,  raised  to  the  second  power.  The  average 
variance  is  the  mean  value  of  the  variance  of  all  previous 
records  in  the  run  ensemble.  This  allows  NA4  to  compare 
the  current  gear  vibration  with  the  baseline  of  the  system 
under  nominal  conditions.  NA4  is  given  by  the  quasi- 
normalized  kurtosis  equation  shown  below; 


NA4*(M)  = 


NE  fri  -  T)^ 

i=l 


(4) 


where 

r  =  residual  signal 

F  =  mean  value  of  residual  sigiial 

N  =  total  number  of  data  points  in  one  time  record 

i  =  data  point  number  in  time  record 

j  =  time  record  number 

M  =  current  time  record  number  in  run  ensemble 


This  provides  NA4  with  enhanced  trraiding  capabilities,  in 
which  the  kurtosis  of  the  current  signal  is  compared  to  the 
variance  of  the  locked  baseline  signal  under  nominal 
conditions. 


NB4  is  another  parameter  similar  to  NA4  that  also  uses 
die  quasi-normalized  kurtosis  given  in  equation  (4).  The 
major  difference  is  that  instead  of  using  a  residual  signal, 
NB4  uses  the  envelope  of  the  signal  banpassed  about  the 
mesh  frequency.  Again,  as  with  NA4*,  NB4*  is  an  enhance¬ 
ment  to  the  NB4  parameter,  in  which  the  value  of  the 
average  variance  is  "locked"  when  the  instantaneous 
variance  exceeds  a  pre-determined  value.  The  equation  for 
NB4*  is  given  below: 


NB4*(M)  = 


N 

NE  (Sj-  s)'* 
i=l 


(5) 


and 

s(t)  =  magnirnde  of  |b(t)  +  i  {H|b(t)]}| 


(6) 


where 

b(t) 

H[b(t)] 

N 


j 

M 


=  time  averaged  signal  bandpassed  filtered 
about  the  meshing  frequency 
=  the  Hilbert  Transform  of  b(t) 

=  total  number  of  data  points  in  one  time 
record 

=  data  point  number  in  time  record 
=  time  record  number 

=  current  time  record  number  in  run  ensemble 


(C)  Joint  Time-Frequency  Technique 

To  examine  the  vibration  signal  in  a  joint  time- 
frequency  domain,  the  Wigner-Ville  method^”^’^^  is  used  in 
this  study.  While  the  FFT  technique  (eq.  1)  can  provide  the 
spectral  contents  of  the  time  signal,  it  cannot  distinguish 
time  phase  change  during  a  con^lete  cycle  of  operation.  In 
other  words,  it  assumes  that  the  time  signals  are  repeatable 
for  each  time  data  acquisition  window  without  considering 
the  effects  of  any  magnitude  and  phase  changes  during  the 
sampling  period.  The  Wigner-Ville  distribution  will  provide 
an  interactive  relationship  between  time  and  frequency 
during  the  period  of  the  time  data  window.  The  WVD 
(Wigner-Ville  Distribution)  can  be  written  as: 


W(t,  f)  =  E  xjt  *  x*jt  - 


(7) 


or  in  a  discrete  form  as 


L 

W^(nT,  f)  =  2T  E  x(nT4  iT)  x*(nT-  iT)  .  (8) 

i=-L 


An  enhancement  to  this  parameter  is  given  by  NA4*,  in 
which  the  value  of  the  averaged  variance  is  "locked"  when 
the  instantaneous  variance  exceeds  a  pre-determined  value.^ 
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where 

W(t,  f)  =  the  Wigner-Ville  distribution  in  both  the 
time  domain  t  and  frequency  domain  f. 
x(t)  =  the  time  signal 

T  =  the  sampling  interval 

L  =  the  length  of  time  data  used  in  the  transform 

To  allow  sampling  at  the  Nyquist  rate  and  eliminate  the 
concentration  of  energy  around  the  frequency  origin  due  to 
the  cross  product  between  negative  and  positive  frequen¬ 
cy,*’^  the  analytic  signal  was  used  in  evaluating  the  WVD. 
The  analytic  signal  s(t)  is  defined  as 

s(t)  =  x(t)  +  jH[x(t)]  (9) 

where 

H[x(t)]  =  the  Hilbert  transform  of  x(t)  defined  by: 

H[x(t)]  =  1  d^  (10) 

jc  t-c 

*  oo  ^ 

However,  an  alternative  approach  can  be  used  to  calculate 
the  analytic  signal  using  the  frequency  domain  definition. 
The  analytic  signal  s(t)  can  be  evaluated  by  calculating  the 
FFT  of  the  time  signal  x(t),  then  setting  the  negative 
frequency  spectrum  to  zero.  The  analytic  signal  can  be 
obtained  by  evaluating  the  inverse  FFT  of  the  spectrum.  To 
simplify  the  computational  effort,  the  WVD  can  be  evaluat¬ 
ed  using  a  standard  FFT  algorithm.  Adopting  the  convention 
that  the  sampling  period  is  normalized  to  unity,  equation  (8) 
can  be  rewritten  as 

L 

W^(n,  f)  =  2  £  x(n*  i)  x*(n-  i)  .  e  (H) 

i=-L 

As  for  the  continuous  time  case,  it  is  necessary  only  to 
evaluate  the  WVD  at  time  zero. 

Hence 

L 

W^(0,  f)  =  2  J]  k(i)  (12) 

i=-L 

where 

k(i)  =  s(i)  s’''  (-i). 

Equation  12  can  be  evaluated  using  the  discrete  FFT 
algorithm. 

In  order  to  avoid  a  repetition  in  the  time  domain  WVD, 
a  weighting  function*^  is  added  to  the  time  data  before  the 


evaluation  process.  Such  a  process  may  decrease  the 
resolution  of  the  distribution,  but  it  will  eliminate  the 
repetition  of  peaks  in  the  time  domain  and  the  interpretation 
of  the  result  will  be  substantially  easier. 

HI.  Description  of  Experimental  Procedure 

The  fatigue  damage  on  the  test  pinion  shown  in  Figures 
1  to  8  was  obtained  using  the  spiral  bevel  gear  fatigue  test 
rig  illustrated  in  figure  9,  at  the  NASA  Lewis  Research 
Center.  The  primary  purpose  of  this  rig  is  to  study  the 
effects  of  gear  tooth  design,  gear  materials,  and  lubrication 
types  on  the  fatigue  strength  of  aircraft  quality  gears.^ 
Because  spiral  bevel  gears  are  used  extensively  in  helicopter 
transmissions  to  transfer  power  between  nonparallel  inter¬ 
secting  shafts,  the  use  of  this  fatigue  rig  for  diagnostic 
studies  is  extremely  practical.  Vibration  data  from  an 
accelerometer  mounted  on  the  pinion  shaft  bearing  housing 
was  captured  using  a  personal  computer  with  an  analog  to 
digital  conversion  board  and  anti-aliasing  filter.  The 
12-tooth  test  pinion,  and  the  36-tooth  gear  have:,  0.5141  in 
pitch,  35  degree  spiral  angle,  1  in.  face  width,  90  degree 
shaft  angle,  and  22.5  degree  pressure  angle.  The  pinion 
transmits  720  hp  at  nominal  speed  of  14,400  rpm.  The  test 
rig  was  started  and  stopped  several  times  for  gear  damage 
inspection.  The  test  was  ended  at  17.79  operational  hours 
when  a  broken  portion  of  a  tooth  was  found  visually  during 
one  of  the  shutdowns. 

rV.  Discussions  of  Results 

A  series  of  pictures  showing  the  deterioration  of  the 
pinion  teeth  at  various  stages  of  the  test  are  illustrated  in 
figures  1  to  8.  In  figure  1  the  initiation  of  a  small  pit  on  one 
of  the  pinion  teeth  during  the  first  shutdown,  at  about  five 
and  a  half  hour  into  the  test  is  shown.  As  the  test  pro¬ 
gressed,  the  rig  was  shut  down  seven  more  times  to  exam¬ 
ine  the  severity  of  the  pitting  and  its  relationship  with  the 
corresponding  vibrations.  Figures  2  to  4  show  the  increase 
of  the  damaged  area  at  the  pinion  tooth  as  the  elapsed  time 
increased  to  6.55,  8.55,  and  10.03  hr  respectively.  Note  that 
in  figure  5,  as  the  el^sed  time  increased  to  12  hr,  the 
danoage  of  the  pinion  tooth  increased  to  75  percent  of  the 
tooth  surface.  At  this  stage,  pitting  also  initiated  on  the 
adjacent  tooth  and  continued  to  grow  as  the  time  increased 
to  14.53  hr,  figure  6.  At  16.16  hr,  the  damage  has  grown  to 
three  adjacent  teeth  as  shown  in  figure  7.  The  test  was 
terminated  when  a  breakage  is  detected  on  one  of  the  three 
heavily  pitted  teeth  at  17.79  hr,  as  shown  in  figure  8. 

Figure  10  depicts  the  running  speed  of  the  test  rig 
during  various  stages  of  the  experiment.  Note  that  there  is 
some  fluctuations  present  in  the  running  speed  after  each 
shutdown,  with  a  magnitude  of  approximately  6  percent 
about  the  nominal  pinion  speed  of  14,400  rpm.  There  is  a 
sharp  change  in  speed  at  approximately  8.75  hr.  These 
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variations  in  speed  create  a  substantial  effect  on  the  vibra¬ 
tion  signal,  which  is  amplified  in  the  NA4  (fig.  13),  NB4 
(fig.  14),  and  the  WVD  (fig.  21)  analysis. 

Figure  11  shows  the  results  of  the  FMO  analysis.  As 
seen  firom  the  figure,  FMO  shows  only  moderate  changes  as 
the  damage  starts  and  progresses.  It  does  not  provide  any 
indication  as  the  damages  extends  to  the  adjacent  teeth, 
resulting  in  pinion  tooth  ftacture.  The  majority  of  the 
variations  in  the  FMO  parameter  are  most  probably  due  to 
the  speed  changes  experienced  during  the  test. 

Results  from  the  FM4  parameter,  as  seen  in  figure  12, 
shows  a  possible  reaction  as  the  pitting  started  to  occur, 
however,  it  does  not  provide  any  coherent  indication  of  the 
severity  of  the  pitting  as  the  damage  increased.  In  addition, 
it  does  not  provide  any  information  to  distinguish  the  pitting 
of  a  single  or  multiple  teeth. 

Results  fi-om  NA4  and  NA4*  are  illustrated  in  fig¬ 
ure  13.  It  is  obvious  firom  the  figure  that  both  NA4  and 
NA4*  provide  a  very  good  indication  of  the  pitting  develop¬ 
ment  on  the  pinion  tooth.  The  magitimde  of  the  parameter 
increases  to  a  nondimensional  value  of  7  after  shutdown  #2 
at  6.55  hr,  and  further  to  a  value  of  17  when  the  pitting 
covers  75  percent  of  the  tooth  surface  at  8.5  hr.  As  expect¬ 
ed,  the  "locked"  denominator  in  NA4*  provides  a  more 
robust  indication  as  the  pitting  progresses.^’^*  Again,  both 
parameters  are  very  settsitive  to  the  speed  variations, 
especially  after  shutdown  #3  and  #6. 

The  NB4  and  NB4*  parameters,  as  shown  in  figure  14, 
show  a  very  similar  trend  to  those  of  NA4  and  NA4*,  with 
a  more  robust  indication  to  the  severity  of  the  damage. 
However,  both  NA4  and  NB4  did  not  provide  any  type  of 
indication  as  the  damage  spread  to  other  teeth,  and  finally 
as  tooth  firacture  occurred. 

Figure  15  shows  a  WVD  (and  corresponding  intensity 
scale)  and  below  the  WVD  the  frequency  spectrum  (firom 
FFT  analysis)  of  a  uruform  sine  wave  signal  shown  at  the 
left  side  of  the  figure.  Note  that  only  the  frequency  com¬ 
ponent  of  the  input  frequency  is  detected  by  both  the 
frequency  domain  analysis  and  the  joint  time-frequency 
distribution.  The  WVD  does  not  exhibit  any  changes  during 
the  1  cycle  (0  to  360  degrees)  rotation  of  the  shaft.  When 
a  short  term  airqtlitude  and  phase  change  is  added  to  the 
system,  as  shown  in  figure  16,  the  frequency  spectrum 
remains  virtually  unchanged.  The  WVD  shows  a  dramatic 
change  of  the  energy  distribution  pattern  at  the  location 
where  the  change  occurs.  The  lighter  shades  of  the  distribu¬ 
tion  display  intficates  a  smaller  vibration  anq)litude,  which 
is  shown  by  the  time  signal  at  the  left  side  of  the  figure. 
Such  an  effect  could  be  possibly  caused  by  a  chipped  or 
cracked  tooth.  Figure  17  shows  the  effects  of  a  short  term 
amplitude  increase  in  the  time  signal  to  simulate  vibrations 
caused  by  gear  tooth  surface  damage.  Note  that  the  frequen¬ 


cy  spectrum  remains  the  same  showing  only  the  componrait 
of  the  exciting  frequency  while  the  WVD  again  provides  a 
good  indication  of  the  anqilitude  increase  by  the  widening 
of  the  shaded  area  to  a  diamond  shape  at  the  corresponding 
"damaged"  tooth  location.  Figure  18  shows  the  effects  of  a 
timP!  decaying  short  tom  amplitude  and  phase  change 
signal.  The  WVD  shows  a  half  diamond  shape  of  shaded 
area,  similar  to  that  of  figure  16,  at  the  location  where  the 
amplitude  and  phase  changes  are  presented.  As  seal  in 
figure  18,  the  ft^uency  spectrum  gives  very  little  indication 
of  the  signal  change. 

Figures  19  to  26  show  the  WVD  and  firequency  spectra 
of  the  spiral  bevel  pinion  vibration  at  various  stages  of 
damage,  corresponding  to  the  photographs  givai  in  figures  1 
to  8.  Figure  19  shows  the  occurrence  of  the  initial  pitting  at 
around  200  degrees  from  the  triggering  point  of  the  gear,  at 
5.5  hr,  of  running  time.  At  6.55  hr,  the  pitting  on  the  tooth 
surface  progressed  to  a  more  noticable  stage,  as  shown  in 
figure  2,  the  WVD  pattern,  figure  20,  begins  to  adopt  those 
of  a  short  term  amplitude  and  phase  change  as  illustrated  in 
figure  18.  Due  to  the  speed  increase  at  8.55  hr,  the  overall 
W^VD  amplitude  increases  substantially  as  shown  in  figure 
21.  At  the  mnning  time  of  10.03  hr,  the  corresponding 
WVD  in  figure  22  shows  the  initiation  of  a  cross  pattern  as 
surface  pitting  in  the  damaged  tooth  becomes  mote  pro¬ 
nounced  which  can  also  be  evidented  from  the  photogr^h 
of  the  pinion  gear  in  figure  4.  This  phenomenon  is  also 
evident  in  the  frequency  spectrum  with  the  existence  of 
sideband  components.  At  12.03  hr,  when  the  pitting  on  the 
pinion  tooth  hax  extended  to  about  75  percent  of  the  tooth 
surfrce  as  shown  in  figure  5,  the  WVD  pattan  exhibits  a 
solid  cross  pattern  extending  over  the  mesh  frcqeuncy  and 
several  of  its  adjacent  sidebands.  The  high  concentration  in 
the  WVD  energy  and  the  initiation  of  a  second  cross  pattern 
at  14.53  hr,  figure  24,  shows  the  advancement  of  the  pitting 
process  on  second  tooth,  as  illustrated  in  figure  6.  This  is 
farther  confirmed  by  the  large  amplitude  of  sideband 
component(above  mesh)  in  the  frequency  spectrum.  At 
16.16  hr,  as  seen  in  figure  25,  the  WVD  pattern  changes, 
showing  more  advanced  damage  pattern  sunilar  to  tiie 
multiples  of  the  decay  of  a  single  short  term  anq)litude 
increase  and  phase  change  demonstrated  in  figure  18.  Such 
phenomenon  is  due  to  the  pitting  of  three  consecutive  teeth 
in  the  pinion  gear  as  shown  in  the  picture  given  in  figure  7. 
The  frequency  spectrum  in  figure  25  shows  a  substantial 
amplitudes  increase  in  the  sideband  components.  The 
discontinuity  of  the  WVD  at  the  mesh  frequency,  shown  in 
figure  26,  similar  to  the  example  shown  in  figure  16  due  to 
the  short  term  phase  change,  is  probably  the  result  of  the 
instantaneous  phase  change  caused  by  the  fractured  tooth,  as 
illustrated  in  figure  8.  The  two  cross  patterns  in  the  WVD 
is  very  distinct  as  the  effects  of  pitting  at  the  teeth  adjacent 
to  the  fracture  tooth  become  more  pronounced.  Note,  also 
that,  as  given  in  figure  26,  the  amplitude  of  the  sideband 
frequencies(above  and  below  the  mesh  frequency)  increase 
substantially. 
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V.  Conclsuions 

Based  on  the  results  of  the  application  of  the  various 

aforementioned  methods,  the  following  conclusions  can  be 

made: 

1)  The  FMO  parameter  shows  only  moderate 
changes  as  the  damage  starts  and  progresses. 
It  also  fails  to  indicate  the  fracture  of  the 
pinion  tooth. 

2)  The  FM4  parameter  shows  a  possible  reaction 
to  the  start  of  the  pitting  process,  however,  no 
coherent  indication  is  provided  for  the  growth 
and  severity  of  the  pit 

3)  The  NA4*  and  NB4*  parameters  show  good 
reactions  to  the  initial  pitting  damage  and 
very  nice  indications  for  the  growth  and 
severity  of  the  pitting  damage.  However  their 
indications  for  the  tooth  fracture  is  somewhat 
unclear. 

4)  The  WVD  provides  vital  information  con¬ 
cerning  both  the  severity  and  the  location  of 
the  pitting  process  in  the  gear  system. 

5)  The  fracture  of  the  gear  tooth  and  its  exact 
location  can  be  pinpointed  using  the  WVD 
technique.  However  a  machine  vibration 
signature  database  is  required  to  interpret  the 
resulting  WVD. 

6)  The  occurrence  and  the  severity  of  gear  tooth 
failure  can  be  reliably  detected  using  a  com¬ 
bination  of  the  time  averaging,  the  frequency 
analysis,  and  the  WVD  techniques. 
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Figure  1 . — Photograph  of  pinion  damage  at  t  =  5.50  hr. 


Figure  2. — Photograph  of  pinion  damage  at  t  =  6.55  hr. 


Figure  3. — Photograph  of  pinion  damage  at  t  =  8.55  hr. 


Figure  4.— Photograph  of  pinion  damage  at  t  =  10.03  hr. 


Figure  5. — Photograph  of  pinion  damage  at  t  =  1 2.03  hr.  Figure  6. — Photograph  of  pinion  damage  at  t  =  1 4.53  hr. 


Figure  7.— Photograph  of  pinion  damage  at  t  =  1 6.1 6  hr.  Figure  8.— Photograph  of  pinion  damage  at  t  =  1 7.79  hr  (end  of  run). 
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-Application  of  WVD  to  the  spiral  bevel  gear  vibration  signal  at  8.55  hr.  Figure  22. — Application  of  WVD  to  the  spiral  bevel  gear  vibration 
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